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Abstract 

By using our novel Grassmann formulation we study the phase 
transition of the spanning-hyperforest model of the /c-uniform com- 
plete hypergraph for any k > 2. The case k = 2 reduces to the 
spanning-forest model on the complete graph. Different k are studied 
at once by using a microcanonical ensemble in which the number of 
hyperforests is fixed. The low-temperature phase is characterized by 
the appearance of a giant hyperforest. The phase transition occurs 
when the number of hyperforests is a fraction (k — l)/k of the total 
number of vertices. The behaviour at criticality is also studied by 
means of the coalescence of two saddle points. As the Grassmann 
formulation exhibits a global supersymmetry we show that the phase 
transition is second order and is associated to supersymmetry breaking 
and we explore the pure thermodynamical phase at low temperature 
by introducing an explicit breaking field. 
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1 Introduction 



The phase transition in a model of spanning forests is particularly 
interesting because only the geometric properties of connection of dif- 
ferent parts are involved and this extremely reduced structure is prob- 
ably at the root of many critical phenomena, within, and even outside, 
natural sciences. 

A possible way to attack this problem on a generic graph by the 
tools of statistical mechanics goes back to the formulation as a Potts 
model [U El |3] in the limit of vanishing number q of states. The Potts 
model on any finite graph G = (V, E), with vertex set V and edge set 
E, is characterized by the coupling v e , for each edge e S E, which is 
related to the inverse temperature (3 and exchange coupling J e by the 
relation v e = e^ Je — 1 . By definition q is a positive integer and the set 
of couplings v = {v e } ee E is of real numbers in the interval [— l,oo]. 
The Fortuin-Kasteleyn representation jH [5] expresses the partition 
function Zc(q\ v) of the Potts model as a sum on all subgraphs HOG 
of monomials in both q and v e 's 

Z G (q;v) := « K(H) U «- (1-1) 

HCG e&E(H) 

where K{H) is the number of connected components of the subgraph 
H. Therefore the model is easily extended to more general values 
of its parameters. In this form it takes the name of random cluster 
model [6] . More generally, it is convenient to introduce the redundant 
description in terms of two global parameters A and p 

Z G (\,p; W ):=Y,* K(H) - K(G) P LiH) II W * (1-2) 

HCG e€E(H) 

where L{H) is the cyclomatic number of the subgraph H. The redun- 
dancy is easily shown by using the Euler relation 

V - K = E - L (1.3) 

and the relations 



q =Xp (1.4) 
v =w p . (1-5) 

Indeed, this form is suitable for taking two different limits when q — ► 
0. In the former limit A — > at p fixed only maximally-connected 
subgraphs will survive and will be weighted by a factor p L ^ H \ In the 
latter one p — » at A and w fixed only spanning forests will survive, 
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weighted by a factor \ K ( H ' as A - ^^ is an overall constant. Remark 
that, when G is planar, maximally-connected subgraphs are in one-to- 
one correspondence with spanning forests by graph duality, and that 
when both A — ► and p — > only spanning trees will survive in each 
connected component of G. 

As shown in [7J |8], the model in the limit p — > 0, that is the 
spanning forest model, admits a representation in terms of fermionic 
fields, which means that the partition function can be written as a 
multiple Berezin integral over anti-commuting variables which belong 
to a Grassmann algebra. Moreover, this representation [HI [9] is power- 
ful enough to describe the model with many body interactions which 
gives rise to hyperforests defined on a hypergraph [10 , which is a 
natural generalization of the concept of graph where the edges can 
connect more than two vertices at once. 

In two dimensions, the critical behaviour of the ferromagnetic 
Potts/random-cluster model is quite well understood, thanks to a 
combination of exact solutions [TT] , Coulomb gas methods [12] , and 
conformal field theory [13] . Information can also be deduced from the 
study of the model on random planar lattices |14|, [TE[ [T6] . Also in 
the q — > limit detailed results are avalaible both for the tree model, 
in particular in connection with the abelian sandpile model [T7|, as 
for spanning forest on a regular lattice [7J [TH] and directly in the 
continuum [19]. Also the model on random planar lattices has been 
considered [2"U] . 

But in more than two dimensions the only quantitative informa- 
tions we have about the spanning-forest model come from numerical 
investigations [21]. Monte Carlo simulations performed at increasing 
dimensionality (d = 3,4,5) show a second-order phase transition. 

Much less results are available for the case of hyperforests. Also 
in two dimensions or in the limit of hypertrees. Even the problem 
of determining whether there exists a spanning hypertree in a given 
/c-uniform hypergraph, is hard, technically NP-complete, for k > 4, 
whereas for k = 3, there exists a polynomial-time algorithm based on 
Lovasz' theory of polymatroid matching [22]. See [23 for a randomized 
polynomial-time algorithm in the case k = 3 whose main ingredients 
is a Pfaffian formula for a polynomial that enumerates spanning hy- 
pertrees with some signs [23] , which is quite similar to our Grassmann 
representation [25] . 

In [26] a phase transition is detected in the random /c-uniform 
hypergraph when a number of hyperedges \E\ = n/k{k — 1) of the 
total number of vertices n = \V\ is chosen uniformly at random. In 
the case of random graphs, that is for k = 2, Erdos and Renyi showed 
in their classical paper [27] that at the transition an abrupt change 
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occurs in the structure of the graph, for low density of edges it consists 
of many small components, while, in the high-density regime a giant 
component occupies a finite fraction of the vertices. Remark that 
their ensemble of subgraphs is the one occurring in the microcanonical 
formulation, at fixed number of edges, of the Potts model at number 
of states q = 1. The connected-component structure of the random 
^-uniform hypergraph has been analyzed in [26] where it has been 
shown that if \E\ < n/k{k — 1) the largest component occupies order 
log n of vertices, for \E\ = n/k(k — 1) it has order n 2 / 3 vertices and for 
\E\ > n/k{k — 1) there is a unique component with order n vertices. 
More detailed information on the behaviour near the phase transition 
when \E\ — > n/k{k — 1) have been recovered in [28, 29 for the case 
of the random graph, but see also [3U1 [3T] , and in [32] for the general 
case of hypergraphs. 

By using the new Grassmann representation, we present here a 
study of the phase transition for the hyperforest model on the k- 
uniform complete hypergraph, for general k, where the case k = 2 
corresponds to spanning forests on the complete graph. The random- 
cluster model [33] on the complete graph has already been developed 
but it cannot be extended to the q — > case, exactly like the mean-field 
solution for the Potts model |34[ [35] . The fermionic representation, 
instead, describes the Potts model directly at q = as it provides 
an exact representation of the partition function of the spanning- 
hyperforest model. 

As usual with models on the complete graph, the statistical weight 
reduces to a function of only one extensive observable, which here is 
quadratic in the Grassmann variables. Under such a condition the 
partition function can be expressed as the integration over a single 
complex variable in a closed contour around the origin [9]. Indeed, 
counting the spanning forests over a complete (hyper-)graph is indeed 
a typical problem of analytical combinatorics. And, exactly like in the 
case of ordinary graph, when the number of connected components 
in the spanning forests is macroscopic, that is a finite fraction of the 
number of vertices, there are two different regimes, which can be well 
understood by means of two different saddle points of a closed contour 
integration over a single complex variable as presented in [36 (but see 
also the probabilistic analysis in [37J). And even the behaviour at the 
critical point can be studied as the coalescence of these two saddle 
points. 

In this paper we shall first review for reader convenience the Grass- 
mann formulation of the spanning-forest model in Sec. [2] and in Sec. [3] 
how it is possible to recover, in the case of the k- uniform complete 
hypergraphs, a representation of the partition function suitable for 



4 



the asymptotic analysis for large number of vertices n. In this same 
Section we shall also present a full discussion of the saddle points in 
the micro-canonical ensemble, that is at fixed number of connected 
components, and of the associated different phases. We shall see that 
the universality class of the transition is independent from k. And we 
will exhibit the relation with the canonical ensemble. In Sec. [4] we will 
provide an interpretation of the transition as the appearance of a giant 
component by introducing a suitable observable which is sensible to 
the size of the different hypertrees in the hyperforest. 

More interestingly, our Grassmann formulation exhibits a global 
continuos supersymmetry, non-linearly realized. We shall show that 
the phase transition is associated to the spontaneous breaking of this 
supersymmetry. By the introduction of an explicit breaking we shall 
be able to investigate the expectation values in the broken pure ther- 
modynamical states. We shall therefore be able to see in Sec.[5]that the 
phase transition is of second order. This seems at variance with the su- 
persymmetric formulation of polymers given by Parisi and Sourlas [38J 
where it appeared to be of zeroth order. 

2 The spanning-forest model 

Given the complete hypergraph K, n = (V, E) with vertex set V = [n] , 
and complete in the fc-hyperedges for all 2 < k < n, so that the 
hyperedge set E is the collection of all A C V with cardinality at 
least 2, let's introduce on each vertex i £ V a pair of anti-commuting 
variables ipi, ^>i\ 

fyiM = {Mj} = {*i>M = Vi,j£V(G) (2.1) 

which generate the Grassmann algebra A[ipi, . . . , tp n , ^i, . . . , ip n ] of di- 
mension 2 2n . 

Then, for each hyperedge A C E, we define the monomial 

TA-=]J^i, (2.2) 

and, for each indeterminate t, the Grassmann element 

ff := t(l - \A\)t a + TA\i - Yl M-^Mm} • ( 2 - 3 ) 

In [8] it has been shown that the generating function of unrooted 
spanning forest on a generic hypergraph admits the following repre- 
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sentations: given a set of edge weights w = {wa}a&e we have 



= f V n $, ^) expltJ2 Mi + E w AfA \ ( 2 - 5 ) 

X>„(^) exp(-W) (2.6) 

where the indeterminate t plays the role of the parameter A we had 
in the random cluster model formulation, T is the set of hyperforests, 
K{F) the number of connected components in the hyperforest F, that 
is the number of hypertrees, 

J V n $,i/>) := J IJ#i^ ( 2 - 7 ) 

is the Berezin integration and we denoted by —TC, as usual in statistical 
mechanics, the exponential weight. 

The fermionic model we introduced above presents a non-linearly 
realized osp(l|2) supersymmetry [8]. Firstly, we have the elements 
of the sp(2) subalgebra, with 

5ipi = -aipi + -ftpi (2.8) 
Sfa = +aA + pA (2.9) 

where a,/3, 7 are bosonic (Grassmann-even) global parameters. Sec- 
ondly, we have the transformations parametrized by fermionic (Grass- 
mann-odd) global parameters e, e: 

8ij} i = r 1 l 2 e{l-t^ i ) (2.10) 
6$i = rV 2 e(l-t$iil>i) (2.11) 

In terms of the differential operators di = d/dipi and di = d/dipi, the 



transformations (2.9) can be represented by the generators 

x Q = y £ j {^idi-i> i d i ) (2.12) 

iev 

X + = J2$i 9 i ( 2 - 13 ) 

ieV 

X- = J2^ B i ( 2 - 14 ) 
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corresponding to the parameters a, /3, 7, respectively, while the trans- 
formations (2.10) and (2.11) can be represented by the generators 



Q+ = t- 1 / 2 Y,(i-t4>iA)d i (2.15) 
Q- = t- 1/2 Y,( 1 ~ t ^)^ ( 2 - 16 ) 

iev 

corresponding to the parameters e, e, respectively. These transforma- 
tions satisfy the commutation/anticommutation relations 

[X Q ,X ± ] = ±2X± [X+,X-] = X (2.17) 

{Q±,Q±} = ±2X± {Q + ,Q-} = X (2.18) 

[X ,Q±] = ±Q± [X±,Q±]=0 [X±,Q T ] = -Q± 

(2.19) 

Note in particular that X± = Q\ and Xq = Q+Q- + 

3 Uniform complete hypergraphs 

The ^-uniform complete hypergraph is the hypergraph whose vertices 
are connected in groups of k in all possible ways or, alternatively, 
whose edge set E is the set of all fc-sets over the vertex set V. In 
our general formulas for the complete hypergraph we must set the 
weights wa = for all the hyperedges A with cardinality different 
from I .A I = k. In the following we will set all the nonzero weights 
to one, so that we shall restrict ourselves to a simple, one-parameter, 
counting problem. In this case the expansion of the partition function 
in series of t 

Z(t) = Y,Z P t p (3.1) 

p 

provides Z p the number of hyperforests, with all hyperedges with car- 
dinality k, composed by p hypertrees. Please remark that, in the 
/c-uniform complete hypergraph, the number p of hyperforests must 
be such that 

n — p , 

5 " *3T < 3 ' 2 > 

must be an integer. Indeed it is the total number of hyperedges in the 
hyperforest. 
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By definition 



1 



Z v = - (UP) t=0 



pi 



V n ($^)W exp<^ 



E U 

yA^A\=k 



(0) 



(3.3) 
(3.4) 



where 



i&V A:\A\=k 

and (-)t=o is the un-normalized expectation value in the ensemble of 
hypertrees. 

The interested reader can find a full comparison of our approach 
with respect to the standard tools of combinatorics in our previous 
paper [9] . We introduce a mean-field variable 



and we observe that 



U = ipip + (1 - k) 



E u 



(0) (ipip) 
— n 



fc-l 



(fc-1)! 



kl 



jfc-2 



(fc-2)! 



(3.6) 

(3.7) 
(3.8) 



where J is the matrix with 1 on all entries, so that (tp, J-0) = Y2i j ^Pi^j- 
The following lemma then applies: 

Lemma 1 ([9j) Let n be the number of vertices, g and h generic 
functions on the Grassmann algebra, then 



V n {i>,^) (WYe hm [1 + Wg{^)) ■ (3.9) 



By this lemma, Z(t) can be written in terms of the sole mean-field 
variable tpif) 



Z{t) = / V n (ip, i/>) exp itU + n 



k-l 



(fc-1)! 



1 



fe-n 



(k-2)l 



. (3.10) 



In order to perform an estimate for the asymptotic value of the integral 
for large n we recall that for an analytic function / 



dt_ no 

2m 



(3.11) 



where the integration contour in the complex plane is around the 
origin. We have the following complex integral representation form 
for the partition function Z(t) 



Z(t) 



nl 



d£ 1 

2~7ri £ 



n+l 



exp < t 



£ + (l-*0 



exp < n 



fc-1 



(fc-1)! 



1 



fc-i 



(k-2)\ 



(3.12) 



Let us first work at fixed number of hypertrees, in a micro-canonical 
ensemble in the physics terminology. 



Expanding (3.12 ) in powers of t we obtain the number Z p of span- 



ning hyperforests on the complete fc-uniform hypergraph which is the 
number of states in the micro-canonical ensamble 



Z„ 



nl 



pi 



d£ 1 

2ri 1 



n+l 



£ + (!-*) 



exp < n 



fc-1 



(fc-1)! 



1 



fc-i 



(fc-2)! 



. (3.13) 



Since we are interested in obtaining Z p in the thermodynamical limit 
n — > oo also for large values of p, we define p = an with fixed a as 
n,p — > oo. Changing the variable of integration to rj = (fc — 11 
we obtain the following integral expression: 



1_ fcT' 



T(an+ 1) 



fc-1 
fc! 



1(a) 



where 



with 



A(n) 

B(rj) 



1(a) := 
1 — fcry 



2vH 



A(?7) e 



nB(>)) 



A; 



fc-1 



+ alog(l — rf) + 



a — 1 
fc-1 



log 77. 



(3.14) 
(3.15) 

(3.16) 
(3.17) 



Please note that the factor fc — 1 coming from the change of variable in 
the integral is exactly compensated by the fact that a full turn around 
the origin in the rj plane is equivalent to fc — 1 turns of the £ variable. 

Precise estimates of integrals of this kind for n — > 00 can be ob- 
tained by the saddle point method (see [36 for a very complete dis- 
cussion of this method). 
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3.1 The saddle point method 

A saddle point of a function B{rj) is a point 770 where -B'(^o) = 0, it 
is said to be a simple saddle point if furthermore B"{t)q) 7^ 0. In this 
case it is easy to see that the equilevel lines divide a neighborhood of 
770 in four regions where HeB(r]) is alternately higher and lower than 
the saddle point value Re -6(770) • We will refer to the two lower regions 
as the valleys. 

Analogously, a multiple saddle point has multiplicity p if all deriva- 
tives up to B^ p \t]q) are equal to zero while -B( p+1 )(7? ) + 0. In this 
case there are p + 1 higher and lower regions. 



When evaluating Cauchy contour integrals of the form (3.15), sad- 
dle points of B{r]) play a central role in the asymptotic estimate for 
large n. The method essentially consists of two basic ingredients: an 
accurate choice of the contour and Laplace's method for the evaluation 
of integrals depending on a large parameter. 

The contour has to be chosen to pass through a point which is 
a global maximum of the integrand along the contour and that a 
neighborhood of which (the central region) dominates the rest of the 
contour (the tails) as n grows. Since an analytic function cannot 
have an isolated maximum, this implies that the contour should pass 
through a saddle point. 

The existence of a contour surrounding the origin and that crosses 
a saddle point along its direction of steepest descent requires that two 
of its valleys are topologically connected and the region connecting 
them surrounds the origin. 

Once we have a contour, we proceed neglecting the tails and ap- 
proximating the functions A(rj) and B{r]) with their Taylor series 
about the chosen saddle point 77*. Then, after having absorbed the 
factor n into a rescaled variable x = (77 — rf)/n 1 ^ p+l " > (where p is the 
multiplicity of the saddle point), we can easily obtain an asymptotic 
expansion of the integral in inverse powers of n. 

We collect here the first few terms of the asymptotic expansion for 
the case of a simple saddle 



nB( V *) 

I 



yj2imB"{rf 



A( V *) + -C(r,*) + ^D(ri*) + o(^ 



(3.18) 



where the terms in the square brackets with half-integer inverse-power 
of n vanish, and of a double saddle 



nB( V *) 

7~ 



1 5(3)^)3 



71 3 



70,4(7/*) + \c( V *) + -D(rj*) + 0(\ 

77,3 n \773 

(3.19) 
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where the terms in the square brackets with powers n~ with 
integer I, vanish. In these formulae C, C, D, and D are rational 
functions of A(r/*), B{r}*) and their derivatives, whose expression is 
reported in the Appendix|Aj together with the value of the constant 70. 



For our integral (3.15 ) in the large n limit the relevant saddle-point 



equation B'(rj) = has two solutions, r] a and 

Va = r r] b = l- a. (3.20) 

If a 7^ a c = {k — 1) jk the two solutions are distinct and correspond 
to simple saddle points. To understand which one is relevant to our 
discussion we need to study the landscape of the function B(r\) beyond 
the neighborhood of the saddles. 

In our specific case, as illustrated in figures from [T] to [3] when 
a < a c among the two saddles only r\ a is accessible, while, if a > a c , 
only rjf, is so. When a = a c the two saddle points coalesce into a 
double saddle point, thus with three valleys, having steepest-descent 
directions e 3 , with k = 0, 1, 2. Of these valleys, the ones with the 
appropriate global topology are those with indices k = 1 and 2. 

As a first result of this discussion, in order to study the asymptotic 
behaviour of Z an will need to distinguish two different phases, and a 
critical point, upon the value of a being below, above or equal to a c . 

We will name the phases with a smaller and a larger number of hy- 
pertrees, respectively, the low temperature and high temperature phase, 
the reason being that, as we shall see, in the low temperature phase 
there is a spontaneous symmetry breaking and the appearance of a 
non-zero residual magnetization. 

3.1.1 Low temperature phase 

In the case a < a c the relevant saddle point is rj a = 1/k. See Fig. [I) 



Since A{r) a ) = 0, we are in the case in which the leading order of (3.18) 
vanishes and the next order has to be considered. The expansion of 
A(rj) and B(n) in a neighborhood of the saddle r] a is as follows: 



A(r] a + u) ~ -k 2 u + k 3 u 2 + 0(u 3 ) (3.21) 

(3.22) 



B(f] a + u) ~ 1 + a ) logfc + a l g ac + fc Qc 01 ^L. 



ka c aj. 



, 9 1 — a a 
k + -3 



y + 0(n 4 ) 



11 



A 




Figure 1: Contour levels for Re5(r/) when a < a c . More precisely, the figure 
shows the case k = 2 and a = \ol c . The two bold contour lines describe the level 
lines of KeB(rj) for the values at the two saddle points (located at the bullets). 
Darker tones denote higher values of HeB(r/). The crosses and the dotted lines 
describe the cut discontinuities due to logarithms in B{rj). The dashed path sur- 
rounding the origin going through one of the saddle points is an example of valid 
integration contour, and the solid straight portion of the path describes an interval 
in which the perturbative approach is valid. 
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A 




Figure 2: Contour levels for KeB(r]) when a > a c . More precisely, the figure 
shows the case k = 2 and a = \a c . Description of notations is as in figure 1 



A 




Figure 3: Contour levels for Re -£>(??) when a = a c . More precisely, the figure 
shows the case k = 2. Description of notations is as in figure [TJ 
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Using formula (3.18) we obtain for (3.13) the following asymptotic 
expression: 



7 ~ 



IV. 



r(an+l) V2vrn 3 p-2)!] n ^r 



Jfe - 1 



-5/2 



n-2 



(an) c 



lev 



2tt 



p-2)!]^ 



1 



ka 
k-l 



(3.23) 

-5/2 

(3.24) 



where in the second line we used the Stirling formula to approxi- 
mate the large factorial nl. In a previous work [9] we already gave an 
asymptotic formula for the number of forests with a given number p 
of connected components. That formula has been obtained keeping p 
fixed while doing the limit n — > oo, in the notation of this paper this 
means taking a infinitesimal. By setting an — > p in (3.23) and using 

-l 

(3.25) 



a 



a 



n 



r(an + l) anT(an) (p — 1)! 

and then taking the limit a — > 0, we can re-obtain the result in [9] by 
using again the Stirling formula to approximate the large factorial n!: 



7. ~ 



n 



n-2 



1 



1 H[ (A; _2)!]S(P-1)! V fc 



fc-1 



p-i 



(3.26) 



3.1.2 High temperature phase 



When a c < a < 1 the relevant saddle point changes into rjf, = 1 — 
a (see Fig. [2]) where the functions A(rjb + u) and B(rji) + u) can be 
approximated at 0{u A ) with 



A(Vb + u) ~ k 
B(rj b + u) ~ 



a — a c 
1 — a 
1 — a 



u 



ti- 



ll 



a, 



(a 
1 



l) 2 (a-1) 3 (a-1)^ 



l--log(l-a) 



A; 



+ 



1 a — a c u 
a a r 1 — a 2 



1 



+ 



+ a log a 
1 



fc a r (1 



a. 



u 

y 



(3.27) 
(3.28) 



The situation is quite analogous to the previous one, with the excep- 
tion that A^rjb) ^ 0, and using formula (3.18) we obtain for (3.13) the 
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following asymptotic expression: 



7 ~ 



n\a an {k-l) 
T{an + 1) 



fc(l-a)(fc-2)! 

1/2 



2irn(l — a) 



ii 



k-l 

(1— a)n 



1 



(3.29) 



27m 



l-Q 

fc-l 



fc(l-a) (jfe-2)! 



(afc-A: + l) 1/2 . (3.30) 



Remark that the saddle point method cannot be applied when a = 1, 
but if we replace 



e fc- 1 



27rniEf r (»}=f + 1 



(3.31) 



J7 rv, 1 



(3.32) 



for a ~ 1 we get 
as we should. 

3.1.3 The critical phase 

When a is exactly a c = (k — l)/k, the saddle points rjt, and rj a coalesce 
into a double saddle point in rj a in which the second derivative vanishes 
along with the first one. The expansion of A{rj) and B(rj) are as in 
pT2l] )- pT22| ) with a = a c : 



A(rj c + u) 
B(rj c + u) 



-k 2 u + k 3 u 2 + 0{u 3 ) 
1 1 , „ k-2, , 

+ TT, — log(fc — 1) + " ~ log k 



k-l k(k-l) 



A; 3 



u" . k 4 (k-2)u 4 



+ 



(k-l) 2 3 (k- l) 3 4 



k-l 



+ 0(u 5 ) 



(3.33) 
(3.34) 



Using (3.19) we obtain the following result 

e^r (fc=i) ^ 3 i/6 r( 2/3)(fc - I) 4 / 3 



2n 



/?,! 



r (V in + 1 ) [(lfe-2)!]*^ 2vrn 2 / 3 
TO) 3 1 /6r(2/3)(fc-l) 4 / 3 



(k-2)\ 



2irn 2 / 3 



(3.35) 



(3.36) 



This formula, for k = 2, can be, in principle, compared with the result 
presented in |36[ Proposition VIII. 11], but unfortunately there is a 
discrepancy in the numerical pre- factor. 
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3.2 The canonical ensemble 



According to the definition (3.14) for the number of forests with 
p = an trees Z an , by evaluating the integral I defined in (3.15) by 
the saddle-point method, when the saddle point rf[a) is simple and 
thus away from the critical point a c , we get the following asymptotic 
expansion for large number of vertices n: 



■/?! 



k- 1 
kl 



a nB(-q*) 



T(an + 1) 
We define the entropy density s(a) as 

s(a) 



n 



1 , Z a . 
-log — 
n n\ 



(3.37) 



(3.38) 



so that we can recover the partition function Z{t) by a Legendre trans- 
formation 



da Z an t ' 



da exp {n[s(o) + a log t] } 



(3.39) 

(3.40) 
(3.41) 



that can be evaluated for large n once more by the saddle-point 
method. Calling a(t) the mean number of trees at given t, we have: 



s'(a{t)) + logt = 0. 



(3.42) 



From (3.37) we see that s(a) still has an a-dependent leading order 
in n 



i \ , , OL — 1 

s\a) ~ —a log n — a log a + a + log 



k\ 



1 



k-1 

that would shift the solution down to 0. By the rescaling 

t = nt 



+ B( V *(a)) (3.43) 



(3.44) 



which is usual in the complete graph, in order to obtain a correct 
thermodynamic scaling, we can reabsorb this factor. The saddle-point 
equation now reads 



s'(a) + log n + log t = 



(3.45) 
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whose solution is 

a= J 7 7" J ' 7 ; cuo) 



fc =± p-2)!]^r t for t< t c 
I-^ti^t for t>l 



where t c = [(k — 2)1] 1 K k x ). And by inversion 

j_ 

'^[{k- 2)l\ fe-i q for a < a c 




for a > av 



(3.47) 



In the ordinary graph case, this means 



[| t < 1 ~ 1 2a a < | 

a = 7 . „ t = I 1 \ 3.48 

4 Size of the hypertrees 

We have shown in the previous sections that the system admits two 
different phases. We want now to characterize these two regimes. Our 
field-theoretical approach provides us a full algebra of observables, as 
polynomials in the Grassmann fields, which we could study systemati- 
cally. However, it is interesting to note that some of these observables 
have a rephrasing in terms of combinatorial properties of the forests 
(cfr. [8]). Furthermore, we are induced by the results of [21] to in- 
vestigate the possibility of a transition of percolative nature, with the 
emergence of a giant component in the typical forest for a given en- 
semble. 

A possibility of this sort is captured by the mean square size of 
the trees in the forest, as the following argument shows at least at a 
heuristic level. If we have all trees with size of order 1 in the large 
n limit, (say, with average a and variance a both of order 1), then 
the sum of the squares of the sizes of the trees in a forest scales as 
(a + a 1 jd)n. If, conversely, in the large n limit one tree occupies a 
finite fraction p of the whole graph, the same sum as above would 
scale as p 2 n 2 + 0(n). 

Furthermore, it turns out that the combinatorial observable above 
has a very simple formulation in the field theory, corresponding to 
the natural susceptibility for the fermionic fields, as we will see in a 
moment. 

Let's start our analysis with the un- normalized expectation 



t ftMO = t J V n {^, ip) i/j^i exp (-H) (4.1) 
= Z(t) (4.2) 



17 



because the insertion of the operator ipiipi simply marks the vertex i 
as a root of a hypertree, and in a spanning forest every vertex can be 
chosen as the root of a hypertree. If we now sum over the index i we 
gain a factor \T\ for each hypertree. Therefore we have: 

t (w)=tj2(^i) = J2 tK{F) J2\ T \H w A= nZ w (4.3) 

iev FeF TeF AeT 

as in each spanning hyperforest the total size of the hypertrees is 
the number of vertices in the graph, that is n. By expanding in the 
parameter t and by taking the p-th. coefficient we get the relation 

1 (4.4) 



Z p (p-1)! 
For the un-normalized two-point function 



(V>i^> = / V n {i)^)^j exp(-W) (4.5) 
we know (see [8]) that 

tm j )= £ t K{F) ^X\w A - (4.6) 

FeF TeF AeT 

i,j connected 

As i and j are connected if they belong to the same hypertree, if we 
sum on both indices i and j we gain a factor |T| 2 for each hypertree 

t j^)> = t 53 = ]r in 2 n w A . (4.7) 

ijeV FeF TeF AeT 

The effect of this observable is to introduce an extra weight for hyper- 
trees in the spanning forests which is the square of its size. 

The average of the square-size of hypertrees in the microcanonical 
ensemble of hyperforests with fixed number p of hypertrees is easily 
obtained from the previous relation by expanding in the parameter t 
and by taking the p-th coefficient, so that 

2 i ((yj,jy;)^~Vo 

The very same method of the preceding Section can be used to 
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evaluate this quantity. Still in a mean-field description, we have: 
(^,^)U p - 1 ) t=0 = (4.9) 



V n {^) (VU^) W^ 1 expert 



k-l 



(k-l)\ 



l-(^JV) 



fc-2 



(fc-2)! 
(4.10) 



n! 



n! 



2vri £ n+1 

- 1 
k\ 



e + (i-fc) 



(fc-l)! 



kip- 1 



exp < n 



fc-i 



(fc-1)! 



1— g 



2vri 



where now 



A(rj) e 



A( v ) 



nB(ri) 



1 



(4.11) 
(4.12) 
(4.13) 

(4.14) 



7/(1 - n)' 

and -B(r/) is the same as before. To evaluate this integral we again use 
the saddle point method. Please note that since the function B{rj) is 
unchanged so are the saddle points. 

Using the general expansion for p = an (3.18) we have 

1 ({ip,Jip)U an - 1 



T(an) 



t=0 A (n*) + ±C(n*) + 0^ 

— = an 



A(n*) + ±C(n*) + 0(±)- 
(4.15) 



Now in the low temperature phase we have A{r] a ) = so in order to 
get the leading term we need C{n*) and as 

A(r ]a ) = -^- and C( Va ) = ^ - % , (4.16) 



k-l 



(4.15) at leading order gives 
{\T\ 2 ) an ~ m 



A(y a ) 
C(Va) 



n 



(a - a c ) 2 ' 



a r — a 



a r 



(4.17) 



so that, as soon as a < a c , a giant hypertree appears in the typical 
forest, which occupies on average a fraction 1 — a/a c of the whole 
graph. In the high temperature instead we have 

1 



a(l — a) 
giving (always at leading order) 

(\T\ 2 ) an ~ ar 



and A(r]b 



k- 



a — a c 
1 — a 



Mvb) _ n 1 
A(rjb) k a — a c 
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(4.18) 



(4.19) 



So that 



lim 



(\T\ 







for a < a c 
for a > a c 



(4.20) 



is an order parameter, but it is represented as the expectation value of 
a non-local operator. We shall see in the next Section how to construct 
a local order parameter. 



5 The symmetry breaking 

In this section we will describe the phase transition in terms of the 
breaking of the global osp(l|2) supersymmetry. According to the gen- 
eral strategy (see for example [39]) let's add an exponential weight 



with an external source h coupled to the variation of the fields (2.10) 



and d2.11| ): 

h 5^(1 -tfafa) = h(n-tW), (5.1) 
The partition function becomes now: 

Z(t,h) = [ Pn(VsV0 e-W^-M"-^). (5_ 2 ) 



We have chosen to add the exponential weight with a minus sign be- 
cause in this way when t is sent to zero with the product h t kept fixed, 
we get, aside from a vanishing trivial factor, the generating function 
of rooted hyperforests. 

More generally, for finite t and h, we have that Z(t, h) can be 
expressed as a sum over spanning hyperforests with a modified weight 

z(t,h)= e n te ~" |T| ( i+/i i r i) ( 5 - 3 ) 

which is always positive, at any n, only for h > 0. 



On the fc-uniform complete hypergraph the partition function (5.2 ) 
is expressed 



Z(t,h)= / V n (tp,tp) exp{tU + htipip} 



exp < — nh + n 



k-l 



(k-l)\ 



1 



W) 



k-11 



(fc-2)! 



(5.4) 



To work in the micro-canonical ensemble we again expand in powers 
oft 

n 

Z(t, h) = Z p {h) t p , (5.5) 

p=0 
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where each term of the above expansion gives the partition function 
at fixed number of components: 



Z p {h) 



pi 



exp < —nh + n 



(k-iy. 



k-l 



2)\ 



-nh 



pi 



{{U + hW) p ) t=0 . 



(5.6) 
(5.7) 



Following the very same steps of the previous section, we can write 
this expression in terms of a complex integral: 



Z p (h) 



nl 



k-l 



with 



r(cm + 1 
I(a,h) :=f£LA(r,)e 



I (a, h) 



nB(n,h) 



where A(rj) is the same as in (3.16) while 
k 



B{n,h) := -h + 



k-l 



rj + a log(l + h — rj) + 



a-1 
k - 1 



(5.8) 
(5.9) 

log??. (5.10) 



I(a, h) can be again evaluated with the same technique as above. 
Let's call r]*{h) the position of the relevant saddle point, which is the 
accessible solution of the saddle point equation 



d_ 
drj 



B( V ,h) 



0. 



V=V*(h) 



(5.11) 



If h > the two solutions are real valued and distinct for every value 
of a and the accessible saddle is simple and turns out to be always 
the one closer to the origin. 

In the following we are going to consider all the functions A, B, C 
and D as evaluated on rf (h) and therefore as functions of the single 
parameter h. 



A(h) = A( v *(h)) 
C(h) = C( V *(h),h) 



B{h) = B(rj*(h),h) (5.12) 
D(h) = D{r/*{h),h) (5.13) 



The asymptotic behaviour of (5.8) is given by the general formula 

e nB(h) 



d3l8l): 



Z an (h) oc 



■s/2imB"[h) 



n 



1 



(5.14) 
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The density of entropy is obtained by taking the logarithm of the 
partition function Z an {h) 



s{a,h) := - log , — . 

The magnetization is then the first derivative of the entropy 
ds 



m(a, h) 



dh 



-1 — a 



-1 — an 



({U + h^) an ) t=0 



an— 1 1 



/i=0 



(5.15) 

(5.16) 
(5.17) 



which is written as the expectation of a local operator, and if we set 
h = in this formula we get 



m(0) = 1 



1 



an— 1\ 



)t=0 



Z an (0) T(an) 







(5.18) 



because of (4.4). In order to evaluate first the limit of large number 
of vertices we use the asymptotic expression for Z an (h) to get 



m(a, h) 



dB{h) 1 1 dB"{h) 



+ 



dh 2nB"(h) dh 



1 8A{h) , i dC(h) 
dh 



n 



<"< (5.19) 



A{h) + lC{h) + 0{^) 



The vanishing of -A(O) in the low temperature phase (a < a c ) has the 
consequence that the two limits n — > oo and h — > do not commute, 
indeed: 



lim lim m(a, /i) 
n— >oo h— >o 



lim lim m(a, /i) 
ft— >o n— >oo 



dB{h) 
dh 



i 



h=0 



h=0 



(7(o) a/i 



h=0 



a, . — a 



> 0. 



(5.20) 



(5.21) 



Remark that the magnetization m vanishes at the critical point lin- 
early and not with critical exponent 1/2 as it is common in mean-field 
theory, the reason being that here the order parameter is not linear 
but quadratic in the fundamental fields. 

In the high temperature phase A(0) ^ and the two limits above 
commute. 



lim lim m(a, h) = lim lim m(a, h) 

n— »oo h— >0 h— >0 ra— >oc 



dB(h) 
dh 



0. (5.22) 



h=0 
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In the study of phase transitions the thermo dynamical limit n — > 
oo has to be taken first. Indeed, the ergodicity is broken in the thermo- 
dynamical limit first and then a residual spontaneous magnetization 
appears even when the external field vanishes. Remark that both the 
free energy and the magnetization vary continuously passing from one 
phase to the other. 

The longitudinal susceptibility xl 



XL(a,h) 



d 2 s(a, h) 



dh 2 



■ a (an — 1) 



(WO 2 (u + hw) 



an—2t 



)t=0 



((U + h^) an ) t=0 



(5.23) 

n [1 — m (a, h)] 2 
(5.24) 



can be obtained from the magnetization: 



XL(a,h) 



dm(a, h) d 2 B(h) 



dh 



dh 2 



+ 



1 



1 



1 1 d 2 B"(h) 1 



2n B"(h) 



n 



dh 2 

8A(h) 
dh 



+ 



d 2 A(h) 
dh 2 



2n B"(h) 2 
+ 



1 d 2 C{h) 
n dh 2 



dB"{h) 
dh 



n 



A(h) + ±C(h) + 0(^ 



+ 



1 dC{h) 
n dh 



A(h)+ l -C(h) + 0(^ 



+ o 



11- 



and taking the two limits in the appropriate order we get 

-l 

I 

-i 



lim lim xl{o-, h) 

h— >Qn^oo 



d 2 B{h) 
dh 2 



u(l-a) ( ac 



h=0 



a 

1-a 



(5.25) 



a < a c 

a > ol c 
(5.26) 



which shows that the susceptibility is discontinuous at the transition, 
with a singularity x( a ) ~ \ a ~ ^cl -1 ) so that the transition is second 
order. Remark that the longitudinal susceptibility appears to be neg- 
ative. This means that in our model of spanning hyperforest there 
are events negatively correlated. It is well known that in the model 
of spanning trees on a finite connected graph the indicator functions 
for the events in which an edge belongs to the tree are negatively 
correlated. This is proven by Feder and Mihail |40j in the wider con- 
text of balanced matroids (and uniform weights). See also [Jl] for a 
purely combinatorial proof of the stronger Raileigh condition, in the 
weighted case. The random cluster model for q > 1 is known to be 
positive associated. When q < 1 negative association is conjectured 
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to hold. For an excellent description of the situation about negative 
association see [H] . 

Still following the analogy with magnetic systems, let us introduce 
the transverse susceptibility 

XT( "' h) " nf{an) ^ 



which, by comparison with (4.15), provides, at h = 

2 ..- l2 



XT (a,0) = -(\T\') an . (5.28) 



n 



In Appendix [B] we prove the identity 

m(a, h) = ^ Xt{ol, h) . (5.29) 

This relation is the bridge between the average square-size of hyper- 
trees and the local order parameter. 

At finite n, when the symmetry-breaking field h is set to zero, we 

get 

m(a,0) = (5.30) 
in agreement with formula (14.41) and 



= -2 XL (a,0) (5.31) 

h=0 



m(a,h) dm(a,h) 
Xr(a, 0) = lim 2 v ' ' = 2 



h^o h dh 



which should be compared with the analogous formula for the O(N)- 
model where it is 

XT (a,0) = (N-l) X L(a,0) (5.32) 

and, in our case, as the symmetry is osp(l|2), N should be set to —1 
as we have one bosonic direction and two fermionic ones which give a 
negative contribution. 

The leading n contribution is 

lxT(a,0) = - XL (a,0) = { n {^y fOT ^ ac (5.33) 

r — — for a> a c 

But, for a < a c , if we first compute the large n limit and afterwards 
send h — ► 0, we know that we get a non-zero magnetization and there- 
fore the transverse susceptibility diverges as 

XT «, 0) ~ 2 - V-^ = t — (5-34) 

h n a c 
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which corresponds to the idea that there are massless excitations, 
Goldstone modes associated to the symmetry breaking. Remark that, 
at finite h, the transverse susceptibility does not increase with n, which 
shows that the average square-size of hypertrees stays finite. 
The longitudinal susceptibility instead 



1 a (1 — a) 
m(a,0) a 2 . 



XL(a, 0) ~ -^t^ V 2 ' ■ (5.35) 



diverges only at a = a c , when the magnetization vanishes. 



6 A symmetric average 

At the breaking of an ordinary symmetry the equilibrium states can 
be written as a convex superposition of pure, clustering, states, which 
cab be obtained, one from the other, by applying the broken symmetry 
transformations. The pure state we have defined in this Section uses 
a breaking field in the only direction we have at disposal where the 
Grassmann components are null. A more general breaking field would 
involve a direction in the superspace to which we are unable to give 
a combinatorial meaning. However, if we take the average in the 
invariant Berezin integral of these fields we give rise to a different, 
non-pure but symmetric, low-temperature state. 
In this Section we shall set t = 1. 

The most general breaking field, with total strength h, but arbi- 
trary direction in the super-space, would give a weight 

n 

fe£[A(l-W j )+^ i +^£] (6.1) 

i=i 

where (A; e, e) is a unit vector in the 1|2 supersphere, i.e. e and e 
are Grassmann coordinates and A is a formal variable satisfying the 
constraint 

A 2 + 2 ee = 1. 
Let us introduce the normalized generalized measure 

dn := d\dede5 (X 2 + 2ee- l) . (6.2) 

A symmetric equilibrium measure can be constructed by considering 



25 



the factor 

F[h;$,il>]:= (6.3) 



J dU exp | - h ^ [A (1 - tpitpi) + eipi + ^e] j 
y dede exp jee - fa ^ ["(1 - ee)(l - + #i + ^e] | 



(6.4) 



(6.5) 

= [l - h 2 (tp, 3if>) + h{n- ijjip)] exp [-h{n- ify)] ■ (6-6) 

where only the last expression is specific to our model, but the previous 
are the appropriate expressions for the model of unrooted spanning 
hyperforests on an arbitrary weighted hypergraph. This function is 
symmetric, for every strength h, as it can be easily checked that 

Q ± F = 0. (6.7) 

If we send h — * this factor is simply 1, but if we first take the 
n —* oo limit and then h — > the expectation value of non-symmetric 
observables can be different. 

The partition function is not changed because of the identity (B.7). 
Indeed 

(F) = (e W [-h(n-^)]) = (l) h (6.8) 

for un-normalized expectation values, because of the relation between 
the transverse susceptibility and the magnetization, equation (B.7), 
which is 

= ( [-h 2 (i>, Ji/j) + h{n- exp [-h{n- ) (6.9) 
= h([-h JV) + (n - MO] ) h (6.10) 

for every h, and therefore also for the derivatives with respect to h. 
But consider for example the magnetization. The insertion of the 
given factor F in the un-normalized expectation provides the relation 

((n - W>))h m : = <( n " ^0 F) = ((n - ^))&+ (6.11) 
+ ([-/ l 2 (Vi,J^) + / l (n-^)] f-Aj exp[-/i(n-^)]> (6.12) 
= 2 ((n - $i>)) h - 2 /» ((V>, J^)) fc = (6.13) 
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Similarly 



mm 



{{i>M)F) 



(6.14) 




(6.15) 



(6.16) 



2{{^,^)) h + {(n-^) 2 ) h 

y j (Aipj + ipjipi + (i - - i>jipj))h 



(6.17) 



(6.18) 
(6.19) 



i 1 - f{i,j})h 



(6.20) 



is the total, not-connected, susceptibility, that is the sum of the lon- 
gitudinal and transverse not-connected ones. And also 



<(n - ^) 2 ) s r = - 2 <$> W))h - <(" - ^) 2 )^ (6.21) 



2 + «n - ^) 2 )r = 2 ((ip, m» + an - n?) h (6.22) 



We have found that in the fe-uniform complete hypergraph with n 
vertices, in the limit of large n, the structure of the hyperforests with 
p hypertrees has an abrupt change when p = a c n with a c = (k — 
l)/k. This change of behaviour is related to the appearance of a 
giant hypertree which covers a finite fraction of all the vertices. As 
the number of hyperedges in the hyperforests with p hypertrees is 
(n — p)/{k — 1), this means that this change occurs when the number 
of hyperedges becomes l/k{k — 1), which is exactly the critical number 
of hyperdges in the phase transition of random hypergraphs at fixed 
number of hyperedges [26] • 

If Z{t) is the generating partition function of hyperforests, where 
the coefficient of t p is the total number of those with p hypertrees, 
in the limit of large n there is a corresponding singularity at t c = 



In our Grassmann formulation this singularity can be described 
as a second-order phase transition associated to the breaking of a 



so that 



as it must occur for a symmetric observable. 



7 Conclusions 



np-2)!]- 1 /^-!). 
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global osp(l|2) supersymmetry which is non-linearly realised. The 
equilibrium state occurring in the broken phase can be studied by the 
introduction of an explicit breaking of the supersymmetry. 



A Saddle-point constants 



Let us use the notation 



X 



in) 



d n 



X(rj) 



v=r) * 



For the simple saddle point we report the combinations C and D in 
terms of functions A and 5 



C = 



D = 



245^ 



12^(1) B(2) 5(3) - 12j 4(2)- B (2) + A ( 35(2) -B(4) - 55 ( 3) 



(A.l) 



11525^ 



385A5j 3) + 144 5^,4(4) 



-210 5 (2) 5 ( 2 3) (4A (1) 5( 3 ) + 3A5 (4) ) 

+21 5 2 2) (±0A {2) Bf 3) + 40^)5(3)5(4) + 5AB ( 2 4) + 8AB (3) 5 (5) ) 
-24 5 ( 3 2) (20^(3)5(3) + 15^1(2)5(4) + 6A (1) 5 (5) + A5 (6) ) . 

(A.2) 

For the double saddle points the necessary combinations are instead 



C = 



AB 



(4) 74 



A 



(i) 



R 4/3 4! 
^(3) 



B 



1/3 
(3) 



71 



(A.3) 



jj _ A {Z) 73 | A {2)B{j) 76 | ^(1)^(5) 76 ^(l)^) 79 



5(3) 3! 



B (3) 



2-4! 1 Bf 3) 7 

?3 



R3 
^(3) 



AB (4) B (5) 79 ^ B (4) 712 ^^(6) 76 



R3 
^(3) 



5f 3) 6! 



4!5! 5* 3) 3!(4!) 3 
The constants 7& are given by 

r I. _ 2±_ 

auu e 3! 

'1 + A; 



2(4!) 2 
(A.4) 



1 . / 1+fe 
Ik ■= sin 2vr — — 

7T \ 3 



(3!)¥ 
3tt 



sin 



2vr 





1 + k 



(A.5) 
(A.6) 
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B Ward identities 

As a result of the underlying symmetry, there are relations among the 
correlation functions, called Ward identities [39]. In this Appendix 
we give a more direct derivation of one of them which simply uses 
integration by parts. 
By definition 

U(0 =Z + (l-kf- (B.l) 



fc! 



so that 



du _ 

~d£ ~ (k - 2)1 



(B.2) 



and therefore the un-normalized expectation value of iptp in presence 
of the symmetry breaking is 

2™ exp rw^y- - ht ^ n -^\^^ exv {tu} ■ 

(B.4) 

Perform now an integration by parts 



in 



d£ 1 



2rri £ n+1 



n 1 



-k-l 



(k-2) 



exp < tU + n 



fc-i 



(fc-1)! 

= n£(t,/i) - ht($,Ji>)) h . 
So that 

nZ(t,h) - t(tpip) = ht((ip,Jip)} 
which expanded in series of t implies that 

n {(U + h^)P) t=0 - (i>i> (U + /i^V) p " 1 )t=o = 



ht(n-£)} (B.5) 
(B.6) 



(B.7) 



P 

or for p = an 



h^MW + h^)^ 1 )^ (b. 



1 — a 



t=o 



{{U + /iVV0 Qn >t=o 



a /i 



an— 1\ 



;t=o 



(B.9) 
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which means that, in the microcanonical ensemble, for every h, we 
have 

m(a,h) = -xr(a,h). (B.10) 
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